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We employ a mean-field approximation to study the Ising model with aperiodic modulation of its 
interactions in one spatial direction. Two different values for the exchange constant, J a and Jb, 
are present, according to the Fibonacci sequence. We calculated the pseudo-critical temperatures 
for finite systems and extrapolate them to the thermodynamic limit. We explicitly obtain the 
exponents f3, 8, and 7 and, from the usual scaling relations for anisotropic models at the upper 
critical dimension (assumed to be 4 for the model we treat), we calculate a, v, u//, rj, and rj//. 
Within the framework of a renormalization-group approach, the Fibonacci sequence is a marginal 
one and we obtain exponents which depend on the ratio r = Jb/Ja, as expected. But the scaling 
relation 7 = f3 (8 — 1) is obeyed for all values of r we studied. We characterize some thermodynamic 
functions as log-periodic functions of their arguments, as expected for aperiodic-modulated models, 
and obtain precise values for the exponents from this characterization. 

PACS numbers: 05.50.+q, 64.60.F-, 64.60.-i 
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I. INTRODUCTION 

Non-uniform systems are interesting and important 
both from the theoretical and experimental points of 
view. Experimentally, there are already several tech- 
niques of surface growth PHI] that let one controls the 
layout of the layers in order to follow, for example, an 
aperiodic sequence. On the other hand, many theoretical 
issues may be raised, concerning the behavior of systems 
with random disorder or aperiodic modulations of their 
interactions; it is the last case that concerns us in this 
work. More specifically, our interest is to calculate the 
critical parameters of the Ising model within a mean-field 
framework and to characterize the log-periodic behavior 
of some thermodynamic quantities. 

The interactions of the model we treat can assume one 
between two different values, and are ordered according 
to the Fibonacci aperiodic sequence. For models that 
have a continuous transition in its uniform version, the 
influence of aperiodic modulations on their critical behav- 
ior is determined by the Harris-Luck criterion Q (which 
seems to hold true for models with first-order transition 
as well @). According to this criterion, the Fibonacci 
sequence is a marginal one; several results show that a 
marginal perturbation leads to a dependence of the crit- 
ical exponents on the ratio between the two different in- 
teractions Using the simplest version of a mean- 
field approximation, we confirm these results and expand 
them to include other critical exponents, in order to test 
scaling relations, and to characterize log-periodic oscilla- 
tions. 

The rest of this work is organized as follows. In the 
next section we present some properties of aperiodic se- 
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quences, define the model we treat, and outline the mean- 
field approximation we use. Our results are shown and 
discussed in Section ITO1 and in Section ITVl we summarize 
our findings. 



II. APERIODIC SEQUENCES AND 
MEAN-FIELD APPROXIMATION 

Aperiodic sequences may be used, for example, to 
model quasicrystals Q: interactions vary according to 
the order embodied in the sequences. These are built 
from substitution rules, in such a way that no subset of 
the sequence is ever repeated. In our case, we define 
an Ising model on a hypercubic lattice, of coordination 
number z, given by the Hamiltonian: 



22i JijSiSj 



(1) 



<hJ> 



such that the sum is over nearest-neighbor pairs on the 
lattice, Jij is the exchange constant between spins Si 
and Sj, which can assume the values J a and Jb in a 
particular spatial direction, according to the respective 
letter in the aperiodic sequence, and Si = ±1, Vi. 

In this work we are particularly interested in the Fi- 
bonacci sequence, which is obtained from the substitution 
rules: 



A -> s(A) = AB, B ->■ s{B) = A. 



(2) 



This means that, from one stage of the construction of 
the aperiodic sequence to the next, all letters A are re- 
placed by AB and all letters B are replaced by A. Start- 
ing with the letter A, the first stages of this sequence are: 
A — > AB — > ABA -> ABAAB ->■ ABAABABA. This 
last finite sequence corresponds to the following sequence 
of interaction constants: J aJbJ aJ aJbJ aJbJ A- In one 
of the spatial dimension of the hypercubic lattice (hori- 
zontal, say) the exchange constants follow this sequence, 



2 



while in the remaining perpendicular hypersurface all in- 
teractions assume the same value, which is the same as 
for the succeeding horizontal bonds. An example of a lat- 
tice constructed this way, in two dimensions, is depicted 
in Fig. CD 

One of the interesting theoretical questions one may 
pose is about the influence of aperiodic modulations on 
the critical behavior of the model, when compared to its 
uniform counterpart. For the case of continuous tran- 
sition on the uniform model, the Harris-Luck criterion 
determines whether the introduction of a given aperiodic 
modulation changes the universality class or not [4| . This 
change is determined by the crossover exponent $, given 
by: 



For the Fibonacci sequence: 



$ = l + d a v(u- 1), 



(3) 



where lo is the exponent describing the behavior of geo- 
metrical fluctuations of the sequence (see below), d a is the 
number of dimensions upon which the aperiodic sequence 
acts (d a = 1 in our case) and v is the correlation-length 
critical exponent of the uniform model. When $ > the 
introduction of the aperiodic sequence changes the crit- 
ical exponents from the values assumed for the uniform 
model (the sequence is said relevant in this case) and 
when $ < the critical behavior of the aperiodic model 
is the same as for the uniform one (irrelevant sequence). 
For $ = 0, the sequence is marginal and previous results 
show that the critical exponents are non-universal: they 
depend on the ratio r = Jb/Ja [6]. In the mean-field 
framework, v — 1/2 for the uniform Ising model and the 
crossover exponent reduces to (see Eq. [3]) : 



* = i(l + w). 



(4) 



Therefore, for ui — — 1 the sequence is marginal, which 
is the case for the Fibonacci sequence, as we will shortly 
see. This quantity and others properties of two-letter 
sequences are obtained from their substitution matrix M. , 
which is defined as: 



,»(») „ S U)^ 



M 
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(5) 



s ( i) 

where n i is the number of letters i that are generated 
by applying the rule s(j). Several features of the se- 
quences are determined by the eigenvalues of A4. The 
greatest eigenvalue (Ai ) determines the rate of growth of 
the total number of letters Af, such that Kf ~ A™, n^$> 1, 
where n is the number of iterations in the construction of 
the sequence. The second greatest eigenvalue (A2) deter- 
mines the wandering exponent uj (Eqs. [3]and[4]) through: 



ln|A 2 | 
InAi 



(6) 



such that the fluctuation in one of the letters, g, is given 
by 1: 



(7) 



M 



(8) 



1. 



and Ai = (1 + y5)/2 and A2 = — A^ 1 , such that lo 
as anticipated. 

Therefore, the aperiodic modulation obtained with the 
Fibonacci sequence, within the mean-field approximation 
applied to the Ising model, is a marginal one, when the 
sequence acts on one of the spatial directions. 

We study the present model (Eq. [T]) within the sim- 
plest mean-field approximation. It may be obtained ei- 
ther from the Bogoliubov inequality [10] with a single- 
spin trial Hamiltonian or, in a less rigorous framework, 
from substituting the magnetization m t for the spin 
Due to the aperiodic modulation, the values of mi vary 
along the direction upon which the aperiodic sequence 
acts (although they are the same for a given hyperplanc 
perpendicular to this direction). The system of equations 
one has to solve is: 



m, = tanh [JQ_im, 
i = 1,..,N 



+ (z — 2)K l m l + Kim i+ i + h] , 

(9) 



where Ki = j3Ji, h = j3H, (3 = l/fc^T, H is a uniform 
magnetic field, ks is Boltzmann constant, T is the tem- 
perature, and iV is the number of hyperplanes on the sys- 
tem (or, equivalently, the size of the aperiodic sequence). 



III. RESULTS 
A. Critical temperatures 

The first task is to obtain the critical temperature T c ; 
our strategy is to calculate pseudo-critical temperatures 
for finite systems and extrapolate the data to the ther- 
modynamic limit. Since the transition is expected to be 
a continuous one and our goal is to calculate T c , we can 
expand Eqs. [5] with H = up to first order on the mag- 
netizations: 



Km = 0, 



(10) 



where 
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and 
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(12) 
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r 


fcsTc/ 'Ja 


0.5 


5.2939768858 


0.7 


5.4801586902 


1.0 


6.000000000000038(64) 


1.3 


6.8300746634 


1.5 


7.4992699398 



TABLE I. Extrapolated critical temperatures for some values 
of the ratio r = Jb / Ja- For r = 1 (uniform model) we obtain, 
within error bars, the exact value, ksTc/J — z — 6. 



the interaction constant J for a given r. More precisely, 
J = Ja(pa + tpb)-, where pa and pb are the fraction of 
letters A and B, respectively, on the infinite aperiodic se- 
quence. These fractions are obtained from the entries of 
the eigenvector corresponding to the greatest eigenvalue 
of the substitution matrix. We notice the agreement is 
quite good; the apparent difference for some regions of 
r comes from the fact that we have few data points and 
have made an interpolation of our data. 



with Ki = (z — 2)Ki — 1 and free boundary condi- 
tions (m^v+i = mg = 0). The interaction parameters 
Ki, K2, Km assume the values Ka or Kb according 
to the respective letter on the Fibonacci sequence. Since 
we do not expect the critical exponents to depend on z, 
we have worked only with z = 6 to obtain the critical 
parameters. 

For temperatures greater than the pseudo-critical one, 
the only solution to this system is rh = 0. So, the ma- 
trix K has an inverse, i.e., det(K) 7^ for this region of 
temperatures. Therefore, coming from above, the first 
temperature such that dct(K) = is the pseudo-critical 
temperature. This procedure is applied to systems with 
different linear sizes L (corresponding to the length of 
the aperiodic sequence, N) and extrapolated to L — > 00. 
One has to be sure that the first temperature such that 
det(K) = is actually obtained, since many tempera- 
tures satisfy this criterion below the first one and they 
tend to accumulate close to the pseudo-critical tempera- 
ture as L increases. 

In order to extrapolate our results to the thermody- 
namic limit, we have used the so-called BST extrapola- 
tion [ll[ in two different ways (see below). The errors of 
our evaluations are obtained as usual for the this method 
of extrapolation [ll| . 

Since we expect log-periodic oscillations on mod- 
els with aperiodic-modulated interactions, the pseudo- 
critical temperatures do not converge monotonically to 
the thermodynamic values: on top of an apparent overall 
convergence, there are oscillations on the values for finite 
L. Therefore, we have also applied the BST procedure 
to every other value of the pseudo-critical temperatures. 
Both procedures lead to the same values in the thermody- 
namic limit. In TableUwe show our results for the critical 
temperatures for some values of the ratio r, extrapolated 
from pseudo-critical temperatures obtained for L up to 
121,393 for r = 0.5 and 1.3 and up to 196,418 for the 
other values of r. Note that we show ten decimal figures 
for r ^ 1, which is certainly enough to obtain precise 
values for the critical exponents. For r = 1 we show all 
figures we are able to obtain, since we can compare it to 
the expected value within the mean-field approximation: 
there is an agreement up to 15 decimal figures. 

In Fig. [2] we compare our values for T c with those ob- 
tained in Ref. 0. The quantity T° is the critical tempera- 
ture for a uniform model with the same mean value J for 



B. Magnetization 

Having calculated the critical temperatures, we can 
now obtain, from the original system of equations (Eq. 
[9]), the magnetization for each plane. The goal is to solve 
this system for rrii for different values of the reduced tem- 
perature t = (T — T c )/T c and of the reduced magnetic 
field h (= (3H). In order to accomplish this we have 
tested three procedures: the first one based on the New- 
ton method [12], the second one uses the secant method 
[12] ] and finally the so-called fixed-point method [l~3| . We 
analyzed the convergence time, for large systems and for 
small values of the reduced temperature, and the accu- 
racy (with respect to known results for small lattices). 
The first method is the less precise, the secant method is 
the most efficient for small values of t, and the fixed-point 
method is the most efficient for large lattices. We have 
chosen the last one, to be able to go to larger systems. 

After a predetermined accuracy is achieved, within the 
fixed-point method, we stop the iterations and calculate 
the mean magnetization as: 

1 L 

m(L) = — rrij. (13) 

i=l 

As discussed elsewhere d,H,[l3] , this quantity may be ex- 
perimentally accessible. We now have to extrapolate the 
values obtained for L —> 00. As expected for aperiodic 
modulated models, oscillations occur as depicted in Fig. 
[21 in order to obtain the value of m = m(L — ¥ 00), we 
have used the extrapolation procedure introduced in Ref. 
8. It simply takes the two last pairs of values for m(L) 
and makes a linear extrapolation with each of them. The 
values mi and mi (see Fig. [3]), obtained for 1/L = 0, 
are then the limits of our estimate for m in the thermo- 
dynamic limit. We then take m = (mi + m-i)j2 and the 
error Am = |mi — m2|/2. From Fig. [3] we clearly see 
that this procedure gives an interval for the magnetiza- 
tion that contains the true value in the thermodynamic 
limit, although it overestimates the error. 

The same procedure was employed to obtain the mag- 
netization for a non-zero magnetic field, which is neces- 
sary to calculate the critical exponents S and 7 (see next 
section). 
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C. Critical exponents 

1. Critical exponent ft 

Our first attempt to estimate the critical exponent (3 
was to fit our data, obtained in the thermodynamic limit, 
as explained in the previous section, to a log-periodic 
function: 



m(t)~(-t)V[log 10 (-t)], 



(14) 



where we assume the following form for the function 
Ppog M (-t)]: 

V [log 10 (-t)] ~ {1 + 5 cos [27rClog 10 (-i) + r0]} . (15) 

Therefore, we obtain for the magnetization: 

m(t) = A(-tf {1 + B cos [2wC\og 10 (-t)+r(p]}, (16) 

where 2ir and r = (v5 + l)/2 are convenient constants 
for the fitting. 

Our results for /3, using Eq. I16[ are shown in Table 
HI] second column. The amplitude of the log-periodic 
term is roughly 5 x 10 ~ 3 for all values of r, except r = 1 
(where this term is not present). Two results are worth 
noting: the exponent for r = 1 (uniform model) is known 
to be 1/2; our result, although near this value is not 
consistent with it. Also, the chi-square per degrees of 
freedom (henceforth referred as x 2 ), is much greater than 
1. This shows that our fitting is not a good one for the 
aperiodic models. 



r 


P 


(«) 


(6) 


(c) 


0.5 


0.56872(4) 


0.5683(2) 


0.56824(4) 


0.7 


0.5439(3) 


0.54558(2) 


0.545553(6) 


1.0 


0.489(1) 


0.49989(4) 


0.49984(2) 


1.3 


0.5270(4) 


0.53033(5) 


0.53041(2) 


1.5 


0.5465(3) 


0.54884(4) 


0.54892(2) 



(b). There is a clear oscillation, as predicted by Eq. [T7] 
the mean value of the fitted curve is the exponent /3. 
Note, however, that for values of log 10 (— t) close to — 1 
the behavior departs from the one predicted. Therefore, 
this interval is not in the scaling region and should not 
be used to study the critical behavior. Our fitting is then 
obtained with the data points in the proper interval. For 
comparison, we show the graph o the log-derivative for 
r = 1, in Fig. [5] no oscillation is present but the devia- 
tion from the expected behavior (in this horizontal 
line) is obtained for —t big enough. Results for (3 with 
this procedure are shown in Table [TTl third column. Al- 
though the value for r = 1 does not include the known 
value for the mean-field approximation, it is closer to 
the expected value than for the previous procedure and 
correct up to the third decimal place. Another improve- 
ment with respect to the previous procedure is that the 
values obtained for \ 2 are orders of magnitude smaller: 
they range from 10~ 3 to 10" 1 . The amplitude of the log- 
periodic term is small, as expected (the maximum value, 
for the values of r we studied, is approximately 1CP 2 , for 
r = 0.5) and increases as we move further away from the 
uniform expected 
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As a final check, we have made fittings using Eq. [TBIbut 
now with a restricted interval of the reduced temperature 
t. We have used the interval in which the log- derivative 
behavior is well described by the data. In Fig. @] (a), 
for example, this interval is —4.9 < log 10 (— t) < —2.5. 
The values so obtained of /3 are shown in Table [III forth 
column: although the result for r = 1 is closer to the ex- 
pected value within the mean-field approximation than 
for the first fitting procedure, it is not better than the sec- 
ond one. Also, x 2 has decreased a great deal, when com- 
pared to the first procedure but it is still orders of magni- 
tude greater than for the log- derivative fitting. Therefore, 
we will take as our results for f3 those in Table [TTl third 
column. 

Finally, we would like to stress the excellent agreement 
between our results for this exponent and those in Ref. 
d (see Fig. |BJ}. 



TABLE II. Magnetization critical exponent, /3, as a function 
of the ratio r, obtained: (a) fitting the data to Eq. 1161 (b) 
using the logarithmic derivative, and (c) fitting the data to 
Eq. 1161 but restricting the interval in log 10 (— t). Numbers in 
parenthesis are uncertainties in the last digit. 

To improve our estimates for /3, we have made another 
procedure, which consists in calculating the so-called log- 
arithmic derivative, namely: 

C(t) ~ d ]°. Sw ^ ~ p + flcoflpfcCtogioH) + ~4>], 

(17) 

where it is assumed that B -C 1 in Eq. [15] This deriva- 
tive is numerically obtained and the data is fitted to the 
previous equation. Examples of the type of behavior we 
obtain are depicted in Fig. @]for r — 0.7 (a) and r = 1.5 



2. Critical exponent 8 

In order to calculate the exponent <5, one has to study 
the dependence of the magnetization on the external uni- 
form magnetic field h. As a log-periodic dependence is 
expected, we also have made all three fitting procedures 
described above for this case. Again the best results are 
obtained for the second one. 

More precisely, we assume the dependence of m on H 
to be (a sgn(_ff ) term is present in the following equation 
but we have omitted it, for clarity): 

m{H) = A\H\%{1 + Scos[27rClog 10 \H\+r<f>]}. (18) 

Therefore, the logarithmic derivative is given by (again, 
taken into account that the amplitude of the log-periodic 
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oscillation is small): 

C(H) ~ = i+Bcos[27rClog 10 (19) 

The typical behavior is depicted in Fig. [7] again log- 
periodic oscillations are present and the critical exponent 
5 is obtained from the previous function (Eq. [T9J . 

The critical exponents are shown in Table IIIII The 
mean-field value for r = 1 is 1/3; our numerical eval- 
uations agrees with this result up to the fourth decimal 
place. For the uniform model, as expected, no oscillation 
is present in the logarithmic derivative. Finally, for the 
values of r quoted in Table IIIII B (see Eq. Q1J]) varies 
from 10~ 3 to 1CP 4 and increases as we move away from 
the uniform model. The values of x 2 ( n °t shown) vary 
from 10~ 3 to 10~ 5 for the aperiodic models and equals 
1CP 9 for the uniform case. These results are evidence of 
good fittings. 



r 


1/5 


0.5 


0.37402(2) 


0.7 


0.358745(3) 


1.0 


0.33328(1) 


1.3 


0.34995(2) 


1.5 


0.360770(5) 



TABLE III. Critical exponent 8 as a function of r, for fittings 
to log-derivative functions (Eq. I19|l . Numbers in parenthesis 
are uncertainties in the last digit. 



3. Critical exponent 7 

We have calculated the susceptibility x(t) using two 
different methods. First, for each reduced temperature t, 
we calculate the magnetization for two different (small) 
magnetic fields and perform a numerical derivative to 
obtain x(t). Alternatively, we can differentiate Eq. [9] 
with respect to H and obtain a system of equations with 
Xi(T), i = 1, N, as the variables. Solving for these, we 
can calculate the susceptibility x(t) = J2i Xi{T)/N . 

For the first method, we used the first two procedures 
quoted in the previous subsections, namely: fitting the 
data to the functions 

X (t) =A\t\-" / {l + Bcos[2irClog w \t\+T<f>]} (20) 

and 

^ = rfi ogl \t\ =-^ + B C0S [^ C lo SlO 1*1 + 0]- 

(21) 

But, contrarily to what happened for the two previous 
critical exponents, it was not possible to identify a clear 
log-periodic oscillation for the log-derivative of x(*)- This 
may be due to the importance of more than one harmonic 



in the behavior of this function [15| ; we could not test this 
hypothesis because our data was not enough to obtain 
one period of the log-periodic oscillation. 

Therefore, for the 7 critical exponent we have only 
obtained results from the fitting to a log-periodic function 
as in the previous equation. These results, although not 
as precise as the ones obtained from the log-derivative 
function, should not be off of the correct values by more 
than 0.6%, according to the comparison made for the 
critical exponents /3 and S. Our results are shown in 
Table ITV1 The mean-field value for the critical exponent 
of the uniform case is 1; our evaluation is 0.1 % off. 



r 


7 


leal = 13(5 - 1) 


A 7 (%) 


0.5 


0.9557(8) 


0.9511(5) 


0.5 


0.7 


0.9812(8) 


0.97522(5) 


0.6 


1.0 


1.0010(3) 


1.0000(2) 


0.1 


1.3 


0.994(2) 


0.9851(2) 


0.9 


1.5 


0.9772(7) 


0.9725(1) 


0.5 



TABLE IV. Susceptibility critical exponent as function of r. 
7 stands for the critical exponent calculated using the fitting 
procedure described in the text, ^y ca i stands for the calculation 
using the equality between exponents 7, /3, and 8, and the 
last column shows the percentage difference between the two 
estimates for 7. 

We have also calculated 7 using the usual scaling rela- 
tion 7 = j3(6 — 1), which still holds true for anisotropic 
models (see next subsection), with /3 and S taken from 
the log- derivative fittings. The comparison is in Table 
IIVI note that the discrepancy is 0.9 % for the worst case, 
which confirms our evaluation that the values would not 
be off by much more than 0.6 %. 

As mentioned earlier, another possible method to ob- 
tain the susceptibility is to perform a field-derivative of 
the system of equations for the magnetization (Eqs. [3]), 
in order to obtain a system of equations for Xi- These 
will be given by the solution of this system, in the same 
manner that we did for the magnetization. The results 
are the same as for the previous method, as expected. 
In particular, we were not able to characterize the log- 
periodic oscillations, either. 

4- Other critical exponents 

We now turn to the calculation of other critical expo- 
nents, using the scaling relation for the free energy for 
anisotropic systems. Due to the presence of the ape- 
riodicity in one dimension, we expect different correla- 
tion lengths in the direction of the aperiodic modulation, 
£// ~ t V/l 1 and along the other directions, £j_ ~ t u , with 
q = v II jv ^ 1 [l^ . Assuming the scaling "ansatz" for a 
system in d dimensions (see Ref. [13, where the scaling 
relation is proposed for two-dimensional models) : 

f s (t,h,L) = b-^- 1+ ^f s (b 1 / v ,by\L/b), (22) 
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where f s is the singular part of the free energy, b is the 
rescaling factor, y^ is a scaling exponent, and L is the 
linear size of the lattice. 

From the above equation, one can show, in the usual 
way, the following relations between critical exponents 



way, 

m 

7 = /3(<5-l); a+2/3+7 



2-v{d-\)-v a . (23) 



Therefore, assuming v = 1/2 (since the aperiodic se- 
quence we study is a marginal one |6|) and d = 4, the 
exponents a and ^ assume the values shown in Table 
IVl Note the good accordance with the mean-field values 
for the uniform model (r = 1) and the expected increase 
of v a and decrease in a when we move away from r = 1. 

Assuming a similar scaling form for the two-point cor- 
relation function T(x,y), where x is the distance along 
the aperiodic direction and y is the distance along the 
remaining (d — 1) directions: 



r(x, y ,t)^t^g(x/\t\-^,y/\t\-"), 



(24) 



one can show that: 



d-2 + ri// = 2/3/v/i; d- 2 + rj = 2/3/u, (25) 

where r(x,0,0) - x d-2+v// an d r(0,y,0) - y d ~ 2+r i. 
Therefore, the exponent rj assumes the usual mean-field 
value, namely rj = 0. The values obtained for the expo- 
nent along the aperiodic direction, rj// , are shown in Table 
W\ assuming d = 4, as before. As expected, the value for 
the uniform model is consistent with the known value for 
the mean-field approximation. However, note that the 
value for r = 0.5 is closer to the uniform results than for 
r = 0.7. Since rj// is close to zero and it is obtained from 
v // , which itself is calculated from scaling relations, one 
expect a higher inaccuracy. 



r 


0.5 


0.7 


1.0 


1.3 


1.5 


a 


-0.0877(9) 


-0.06638(9) 


0.0002(2) 


-0.0458(3) 


-0.0701(2) 


v /l 


0.5877(9) 


0.56638(9) 


0.4998(2) 


0.5458(3) 


0.5701(2) 


■nil 


-0.066(4) 


-0.0734(4) 


0.0004(9) 


-0.057(2) 


-0.0746(8) 



TABLE V. Critical exponents calculated from scaling rela- 
tions for anisotropic systems. 



IV. CONCLUSIONS 



We employ a mean-field approximation to treat an 
Ising model with aperiodic modulation in one spatial di- 
rection. The particular aperiodic sequence we use is a 
marginal one, in the context of the Harris-Luck crite- 
rion. We calculate many equilibrium critical exponents, 
including v n and rj//, assuming d — 4 to be the upper 
critical dimension of the model and a particular scal- 
ing form for the singular part of the free energy per site 
and for the two-point correlation function, suitable for 
anisotropic models. As expected, the exponents (with 
the exception of v and 77) depend on the ratio r = Jb/Ja 
but obey the usual scaling relations for anisotropic mod- 
els, whenever possible to test these relations. Our results 
are in accordance with the known values for the mean- 
field procedure (uniform model, r = 1) or with previous 
results for the exponent (3 and critical temperatures [6|. 
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FIG. 2. Comparison between our values for T c (traced line) 
and those obtained in Ref. [6j (continuous line). In the text 
we define the quantity T°. 



FIG. 3. Typical behavior for the magnetization, for fixed 
reduced temperature t or fixed reduced magnetic field h, as 
function of the linear size of the lattice, L. Note the oscilla- 
tory convergence to the thermodynamic limit, as expected for 
aperiodic modulated models. 



FIG. 1. Example of a lattice with an aperiodic modulation 
given by the Fibonacci sequence. In the horizontal direction 
the exchange interactions follow this sequence while in the 
vertical "planes" they assume the same value, equal to the one 
in the following horizontal bonds. Traced (full) lines represent 
J a {Jb) interactions. 



FIG. 4. Logarithmic derivative of the magnetization for r = 
0.7 (a) and r = 1.5 (6). The continuous lines are fittings using 
Eq. [17] while the points are our numerical data. 

FIG. 5. Logarithmic derivative of the magnetization for r = 1. 
The continuous lines are fittings using Eq. [I7]while the points 
are our numerical data. Error bars are approximately the 
same size of the points. 
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FIG. 6. Exponent /3 as a function of r from the data obtained 
in this work (continuous line) and from the results from Ref. 
6] (traced line). 



FIG. 7. Field logarithmic derivative of the magnetization, for 
r — 0.7 (a) and r = 1.5 (b). The continuous lines are fittings 
using Eq. [19] while the points are our numerical data. Error 
bars are approximately the same size of the points 
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